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1. Introduction 



One of the beautiful aspects of the matrix-model formulation of c < 1 string theory is 
that it gives a natural and mathematically precise formulation of the partition function of 
strings moving in different backgrounds. This result began with Kazakov's fundamental 
discovery of the appearance of matter fields in the one-matrix model and culminated in 
the discovery of the generalized KdV flow equations and the associated Wn constraints in 
the c < 1 matrix models coupled to gravity f^-§| . Recently these results have been further 
deepened through the use of a Kontsevich matrix model representation for the tau functions 
relevant to these flows 0, see also ||,|). Analogous results in the c = 1 model have been 
strangely absent, and this paper is a first step in an attempt to change that situation. 
Using recently developed techniques for calculating tachyon correlators in the c = 1 model 
we derive a simple and compact expression (equation (3.10)) for the generating functional 
of tachyon correlators, or equivalently the string partition function in an arbitrary tachyon 
background, valid to all orders in string perturbation theory. In Euclidean space this 
quantity can be interpreted as the partition function of a nonlinear sigma model as a 
function of an infinite set of coupling constants for a set of marginal operators. 

Upon appropriate analytic continuation to Minkowski space the partition function may be 
interpreted as the string S'-matrix in a coherent state basis. 

One immediate consequence of our result (3.10) is that the partition function is natu- 
rally represented as a tau function of the Toda hierarchy. From this result we obtain 
flow equations (equation 4.10) when the c = 1 coordinate X is compactified at the self-dual 
radius. Moreover, this expression can be used to derive a Kontsevich-Penner representa- 
tion of the partition function as a matrix integral, as described in section five below. In 
section six we discuss how time-independent changes in the matrix model background fit 
into our formalism, and in section seven we discuss some open problems and the relation 
of this work to other recent papers on c = 1 and . 

2. Defining correlation functions 

In a particular background, string propagation in a two-dimensional spacetime is 
described on the string worldsheet by the conformal field theory of a massless scalar X 
coupled to a c = 25 Liouville theory cf> with worldsheet action (excluding ghosts) 




dXdX + d(j)d(j) + V2R {2 U + ne 



V2<t> 
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Via its dual interpretation as the conformal gauge action for the coupling of X to two- 
dimensional gravity, (|2.1|) is expressible as the continuum limit of a sum over discretized 



surfaces. The discrete sum, as is by now well known, is generated by a matrix integral. In 
the double scaling limit which leads to the continuum theory this is in turn equivalent to 
a theory of free nonrelativistic fermions with action 

S=y_ oB W^_ + s? _v(A))*. (2.2) 

The potential V(X) in ( |2.2| ) is required to approach — ^A 2 for large A in order to reproduce 
the topological expansion of string theoryi 

The theory contains one field theoretic degree of freedom, the massless "tachyon". 
Tachyon correlators are calculated in the theory (|2.1|) by the insertion of vertex operators 



q ~n-\a\)L ' ( ] 

In this section we show how these correlators are calculated in the double scaled matrix 
model ( |2.2| ). The presentation is a modification of the original derivation in PJ, which 
stressed the spacetime interpretation via collective field theory, in that here we emphasize 
the relation to macroscopic loop amplitudes. 

We recall that the objects calculated in a sum over continuous geometries on two- 
surfaces with boundary are "macroscopic loop amplitudes" defined by fixing the boundary 
values of the two-metric e^ 2 ^ so that the bounding circles C have lengths £ = § c e^/^ 
P,|TT|-|T6|| . In |n],|16] tachyon correlators were defined as the coefficients of nonanalytic 
powers of i in the small-£ expansion of macroscopic loop amplitudes. In particular, in this 
limit the macroscopic loop may be written as a sum of local operators 

W in (£,p) = -r v -^-n-\v\!H^{2^x£) -J2B r , p ^^^ r / 2 I r (2^I£) (2.4) 
sinir \p\ 11 *— ' r z — rr 

1 1 r=l 



2 The behavior of V(A) for negative A is irrelevant to all orders of perturbation theory in 
Indeed, the results of this paper should be interpreted in this perturbative sense. Many results 
are true in the nonperturbative context and we will indicate this in the appropriate places. Where 
we mention nonperturbative results we will refer to potentials which grow sufficiently rapidly for 



large negative A. In [10] these were termed "type I" models. 
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where B r ^ p are redundant operators for p ^ 7L. We may thus extract tachyon correlators 
from macroscopic loop amplitudes as 

n n n 

(n w (^»)>=n r (-i»i)^ l (<n r «> + °w)) + anai y ticin ^ 

i=l i=l i=l 

The matrix model formulation of the theory leads to a simple computation of the appro- 
priate limits of loop amplitudes. In the matrix model the macroscopic loop is related by a 
Laplace transform to the eigenvalue density p(A, x) = i/j*ijj(X, x): 



W{£,x) 



/•OO 

/ e- ex p(\,x)d\ 
Jo 

/oo 
e iqx W{£, x)dx . 
-oo 



(2.6) 



Defining 



we recover 



/•OO 

W(z,x)= e~ zl W{£,x)d£ (2.7) 
Jo 

p(X,x) = --Disc(W(z,x))\ z __ x . (2.8) 

TV 

Inserting ( |2.5| ) we find now 

(f[ p{K qi)) = fl (<fl + °(K 2 )) (2-9) 

i=l i=l i=l 

for generic qi. In the next section we show how this leads to a simple formula for the 
correlation functions. 



3. Calculating correlators 

3.1. Graphical rules for tachyon correlators 

The calculation of tachyon correlation functions now reduces to the study of asymp- 
totics of correlation functions of the eigenvalue density. An explicit formula for these was 
found in [IBJ; we will use a generalization of this to a compactification radius (3 for the 
Euclidean X field. Momenta for the tachyon field are therefore always of the form q = n/[3 
for n G 2. The required modification follows by interpreting the compactified Euclidean 



time as a finite temperature for the free fermion system |T7[]. Allowed fermionic momenta 
are of the form p m = (m+ h)/P with m E 7L. 



The eigenvalue correlator can be written as: 



m= — oo aE'En k=l 



(3.1) 



where E n is the set of permutations of n objects, Q£ = p m + Qa(i) + • • • + q a (k)- The A 
dependence of the correlator is determined by the function I(q, Ai, A2) = (I(—q, Ai, A2))* = 



(A 



1\W- 



-fi-iq 



IA2) where H is the one-body Hamiltonian for the free fermion system. This 



has large A asymptotics given by a direct and a reflected contribution []T0 



I(q, Ai, A 2 



a^+oo — 1 



vAiAl 



D[q,\ u \ 2 ] +R[q,X 1 ,X 2 ] 



D[q,X ll X 2 ] = exp 



i(A* + ig)|log(Ai/A 2 )|- 7 |A?-Ai| 



(3.2) 



R[q,X 1 ,X 2 ] = R q exp 



i(fi + iq) log(AiA 2 ) - 7 ( A i + A l) 



for q > 0. The function R q is the reflection coefficient for the nonrelativistic free fermions 
in the double scaled potential V. For the "standard" case V = —\ A 2 with an infinite wall 
at A = it is given by 



Rn 



1 + ie-^^+iQ) IT{\ -ifj, + q) 
- ie-K^+iq) V r(i +ifi- q ) 



(3.3) 



Inserting ( |3.2|) in the expression ( |3.1| ) leads to a sum of terms. The calculation of tachyon 
correlators requires the extraction of those terms in the sum with the correct asymptotic 
dependence on A^. For each permutation a, at most a finite number of terms in the sum 
over the loop momentum p m contribute to the result. 



A graphical procedure for performing this extraction was developed in |10[ and used 
to derive an explicit expression for arbitrary tachyon correlators. We divide the tachyon 
insertions into "incoming" (q < 0) and 'outgoing" (q > 0) particles. As in a Feynman 
diagram there is a vertex in the (x, A) half-space corresponding to each operator p(x, A). 
While the final result will of course be independent of the order in which the Ai are 
increased to infinity, in intermediate steps we will choose some order and locate the vertices 
accordingly. Points are connected by line segments, representing the integral /, to form a 
one-loop graph. Since the expression for / in ( |3.2|) has two terms we have both direct and 
reflected propagators as in fig. |l|. Each line segment carries a momentum and an arrow. 
Note that in fig. [1] the reflected propagator, which we call simply a "bounce," is composed 



of two segments with opposite arrows and momenta. These line segments are joined to 
form a one-loop graph according to the following rules: 

RH1. Lines with positive (negative) momenta slope upwards to the right (left). 
RH2. At any vertex arrows are conserved and momentum is conserved as time 
flows upwards. In particular momentum q± is inserted at the vertex in fig. [| 
RH3. Outgoing vertices at (x ou t, X ou t) all have later times than incoming vertices 

Diagrams drawn according to these rules correspond to possible physical processes 
in real time and were hence termed "real histories". The connected tachyon correlation 
function is found by summing the terms in (|3.1| ) corresponding to all real histories, and 
reads schematically 

n 

di^)=H) n E ± E n rq(- r qt- ( 3 - 4 ) 

i=l RH m bounces 

The graphical rules allow one to convert ( |3.4| ) into an explicit formula for the amplitude 
||1C|| . In the next subsection we will show that this result may be written quite simply 
in terms of free fermionic fields, representing a fermionized version of the free relativistic 
bosonic field which describes the asymptotic behavior of the tachyon. 



3.2. Free Energy in terms of free oscillators 



One of the central results of [T(| is that the graphical rules described above are equiv- 
alent to the composition of three transformations on the scattering states: fermionization, 
free fermion scattering, and bosonization: if^b°Sffoi b ^f as in fig. The various real his- 
tories correspond to the possible contractions among the incoming and outgoing fermions, 
and the fermion scattering matrix describes a simple one-body process, given essentially 
by the phase shift in the nonrelativistic problem. It should be noted that this does not 
imply the (false) statement that bosonization is exact for the nonrelativistic fermion prob- 
lem. Rather, it is a statement about the asymptotics of certain correlators in the theory 
for a particular class of potentials. Here we will rewrite the tachyon amplitude using this 
formulation as a matrix element of a certain operator in the conformal field theory of a 
free Weyl fermion. 



It is convenient to define rescaled tachyon vertex operators V q = fi 1 ^ q ^ 2 T q and two 
free scalar fields d(fi in / out = J2 n a l n^ out z~ n_1 , such that 



3 = 1 



1=1 
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(3.5) 



where |0) is the standard £X(2,1R) invariant vacuum. The equivalence of the graphical 
rules to bosonization then implies that while the relation between the two bosonic fields 
is complicated and nonlinear it may in fact be expressed as a simple linear transformation 
in the fermionized version. Thus we write dtp = ip(z)tp(z) where ip, ip are Weyl fermions of 
weight | with expansions 

*w = E ^ + ^" m_1 



1Z 



-m— 1 



m6Z 



{VV, IPs} = <W S , • 



Now the result of |T(| states that ( |3.5| ) is equivalent to 



(3.6) 



.out 



— (rn+2) —( m +2> 



.out 



7p m i =R* r i 
Unitarity of the tachyon S'-matrix is equivalent to the identity 



(3.7) 



R q R*_ q — 1 



(3.8) 



on the reflection factors.! Using this, we can rewrite ( |3.7| ) as a unitary transformation 



^ n {z) = S^ out (z)S- 1 
tp tn (z) = StP° ut (z)S- 1 



S =: exp 



•m£Z 



(3.9) 



3 This holds to all orders in perturbation theory for any of the potentials we consider. The 
question of its nonperturbative validity was discussed in p0| . Essentially, this requires that V(X) 
grow sufficiently rapidly for large negative A. 
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Thus we may write the full generating functional for connected Green's functions in terms 
of a single free boson with modes a n : 



i 2 T= ( e S™>i*" Kl/ ' 3+ ^n>i f ™ V '-™ / ' 3 ) ( 



-^(0|e M2 -™>i n n Se M ^->i " "|0) c . 



With this definition T has a genus expansion = JF + A-Fi + • • •• This formula is an 
enormous simplification over previous expressions for c = 1 amplitudes. The generating 
function for all amplitudes is 

Z = e^ T . (3.11) 



4. VFi-i-oo constraints 



In correlation functions of tachyons with integer (Euclidean) momentum, the bounce 
factors R q of ( |373l) simplify due to the following identity 



(4.1) 



This is valid to all orders in perturbation theory. (In the "standard" potential it also 
holds nonperturbatively if n G 2Z.) Note that equation ( |2.9|) holds for generic momenta; 
the results for integer momenta are defined by continuity. Working with the generating 
functional of all amplitudes (|3.11| ) we have: 



i(3 dZ 



ifj, 



/x dt n 
J w n J z ^ Pm 



E„> 1 t » a »( (5a _ n>s -i) > 5 e ^E„> 1 t 



71 *-* — n 



10} 



w , 



n//3-p, 



,.(-) 



(0|e 



(4.2) 



ijj(z)ip(w)Se 



% fjL 



E __ -. ^ ri <^ — n. 
n > 1 
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At the self-dual radius /3 = 1, where all tachyon momenta are integral, we may simplify 
the sum on m using ( fl.l| ) 



d 



,w . 



;-z^-^ + z-) n £(^ 



(4.3) 



the latter sum acting like a delta function. Now integrate by parts and use the identity 

(4.4) 
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(-ill - z^-) n = (-zyz-^i^-Yz^ 



It is convenient to bosonize ip(z) = e^( z \ip(z) = e and shift the zero mode: 



4>(z) = (p(z) + ifilogz . (4.5) 

Taking the operator product of the two exponentials in 0, and using the delta function 
and charge conservation we find the operator: 

1 dw(i^)- n -^— : e~^ w) dl +1 e^ w) : . (4.6) 
/ n + 1 

Now go to the coherent state basis in the i n 's, and redefine the scalar field by a factor of 
i\i to obtain the final result: 

- Z' 1 <b dw^- : e-wM^n+i^M . z ^ ^ 



dt n J n + 1 

where 

d<p(w) = - + V ntnw n ~ l -^—w~ n ~ ^ . (4. 
w ^ a 2 ^ ot n 

n>0 r n>0 n 



The genus zero result of []T8| is easily obtained from this as the leading term at large \i. 
(Note that this was obtained at (3 = oo but genus zero correlators are independent of (3 
Q) 

The operators p( n )(z) =: e - ^ 2 * 1 •Q n e 4>{z) . anc j ^gjj. derivatives generate the algebra 



V^i-i-oo [jT9| . The standard generators are related to these by 



^H^)=E ,i:l£:'';! ^ ( "-' ) w. «.») 



( _ o (n - (2n - 2), 



The rescaling of the scalar field required to obtain ( |4.8| ) is simply a change of basis effected 
by the operator : e 10 ^ 1 ^ 71 "^ : where tt^ is the momentum conjugate to <fi. Inserting this we 
can rewrite (14.71) as 



dZ 



wt +1) z 



di 



(4.10) 

W ( _ n n +1) = & dzW {n+1) (z) . 

Generalizations of ( [4.10|) to other radii (5 ^ 1 follow from (|4.2|) . 

In the conclusions we comment briefly on the relation of this result to the many other 
occurrences of Woo in this subject. 



5. Tau-functions and the Kontsevich-Penner matrix integral 

In this section we will point out that the above reformulation of the generating func- 
tional of the c = 1 string represents mathematically a r-function of the Toda Lattice 
hierarchy. The Toda Lattice naturally contains the KP and KdV hierarchies, and thus 
the c = 1 results are closely related to the expressions obtained for c < 1. We will also 
show how to rewrite the partition function (at the self-dual radius) as a matrix integral, 



generalizing expressions previously considered by Kontsevich and Penner [£0 



5.1. Grassmannians and tau-functions 

Let us first briefly explain the notion of a tau-function and its relation with the 



universal Grassmannian. For more details see e.g. |y] and ||22|1 . We will focus here on the 
relation with conformal field theory instead of the Lax pair formulation. 

Consider a two-dimensional free chiral scalar field ip(z), with the usual mode expansion 

di f (z) = J2^nZ- n - 1 . (5.1) 

n 

The reader is encouraged to think about this scalar field as the target space tachyon field 
at spatial infinity with a periodic Euclidean time coordinate. We have a Hilbert space 
ri built on the vacuum |0), and as in the case of a harmonic oscillator one can consider 
coherent states, 

oo 

\t) = exp ^ ttnOl-n (5.2) 
n=l 

and their Hermitian conjugates 

oo 

(t\ = (0| exp ^ - il nOt n (5.3) 



n=l 



(The parameters t n are considered to be real here.) Now to any state \W) in the Hilbert 
space TC we can associate a coherent state wavefunction Twit) by considering the inner 
product 

r w (t) = (t\W) . (5.4) 

This function is a tau-function of the KP hierarchy if and only if the state \W) lies in the 
so-called Grassmannian. 
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To explain the concept of the Grassmannian we have to turn to the alternative de- 
scription of this chiral conformal field theory in terms of chiral Weyl fermions i^>(z), ip(z) 
by means of the well-known bosonization formulas^ 

idif = V = e itp , $ = e~ i<p . (5.5) 

Loosely speaking, the Grassmannian can be defined as the collection of all fermionic Bogo- 
liobov transforms of the vacuum |0). That is, the state \ W) belongs to the Grassmannian if 
it is annihilated by particular linear combinations of the fermionic creation and annihilation 
operators. 

oo 

^ nH -J2 A nm^- m+ ^\W)=0, n>0, (5.6) 

TO=1 

or equivalent ly, 

\W)=S-\0), S = expJ2 A ^-™-^-mH ■ ( 5 - 7 ) 

n,m 

Note that the operator S can be considered as an element of the infinite-dimensional linear 
group, S e GL(oo,€). 

By replacing the vacuum |0) by the state \ W), we simply made another decomposition 
into positive and negative energy states, and filled these new negative energy states. The 
positive energy wave-functions are no longer of the form z n (n > 0) but are now given by 
the functions 

oo 

V n (z) =Z n ~J2 AnmZ^ . (5.8) 
m=l 

If one prefers the language of semi-infinite differential forms, we have a formula of the form 

\W) =v Q Av x A... (5.9) 

In the above fashion one generates solutions to the KP hierarchy. This construction can 
be extended to give a tau-function for the Toda Lattice hierarchy by considering a second 



set of times t^, as discussed in detail in ||23| . In terms of our conformal field theory, the 
Toda tau-function is simply obtained as 



T 



(t,t) = (t\S\f) , (5.10) 



4 Since we do not wish to flaunt tradition we change conventions for bosonization in this section 
relative to the previous sections. 
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with \t) and (t\ the coherent states ( |5.2| ) and ( |5.3|) and S a general GL(oo,C) element, i.e. 
an exponentiated fermion bilinear of type ( |5.7|) . The S- matrix ( |3.9|) of the c = 1 string is 
definitely of this form, which allows us to conclude that the string partition function can 
indeed (after a rescaling t n — > fi ■ t n ) be identified as a Toda Lattice tau-function. 

Instead of taking the inner product of the state \W) with a coherent bosonic state, 
one can also consider fermionic iV-point functions (see e.g. 0). In fact, one finds in this 
way a simple expression in terms of an iV x N determinant of the wave-functions 



{N\1>{z 1 )...1>(z N )\W) = detv j - 1 (z i ) . (5.11) 

Using the bosonization formulas, one recognizes this correlation function as a special co- 
herent state where the parameters t n are given by 



N 

t 



N 1 

£-*r n - (5-i2) 



n 

i=l 

With this choice of parameterization, and after taking into account a normal ordering 
contribution, the tau-function can be written as 

rit) = , (5.13) 

with A(z) the Vandermonde determinant A(z) = det^ - = rii>j(- 2 i ~ z j)- We apply this 
result in the next section. 

5.2. Kontsevich integrals and the c < 1 models. 

Since our expression for the c = 1 partition function is very analogous to the result 
found for the c < 1 string theories, we will briefly summarize the latter (see [|2|FJ,Q.) Recall 
that the c < 1 matrix models naturally give rise to a universal set of observables O n 
(n = 1, 2, . . .) whose correlation functions 

(O ni ---O ns ) g (5.14) 

at a specific genus g are unambiguously determined. The generating functional r(t) of 
these correlators has an asymptotic expansion in the string coupling constant A 

oo oo 

logr(t) = ^A 2 - 29 (exp^t n O n ) s . (5.15) 

g=0 n=l 
11 



The techniques of the double-scaled matrix models leads to two important results. First, 
the partition function r(t) is a tau- function of the KP hierarchy, that is, it can be written 

as 

r(t) = (t\W) = (t\S\0) , (5.16) 

for some state \W) and matrix S G GL(oo,C). Secondly, all minimal models of type 
(p, q) with fixed p belong to one KP orbit. More precisely, relative to a convenient choice 
of origin, the (p, q) model is obtained at the value tk = Sk,p+q- Furthermore, the KP 
hierarchy reduces to the p th KdV hierarchy, which implies that all correlation functions of 
the operators O n with n = (mod p) vanish. 

The state \ W) corresponding to this orbit is most simply described at the (p, 1) point, 
where a description in terms of topological field theory can be given. For its basis one can 
take the wave-functions 



v n (z) = J^ e ^ P+1 /" ■ Tdyy-. e^-^A , ( 5 .17) 



2ttA 

where the normalization is chosen such that we have the appropriate asymptotic expansion 

v n (z) =z n (l + 0(z- 1 )) . (5.18) 

Since the wave-functions are moments in a Fourier transform, the fermionic formula ( |5.13| ) 
can be explicitly evaluated, and gives rise to a so-called Kontsevich integral [|7| (see also 

]) 

r(t) = c- f DY ■ e «X^-^)A . (5 .i 9 ) 



Here Y and Z are both N x N Hermitian matrices, and the parameterization of the KP 
times tk in terms of the matrix Z is 

A. 



t k = -^TrZ~ k . (5.20) 

This result can be generalized to the 'generalized Kontsevich model' || which features an 
arbitrary potential V(z) 



r{t) = c(Z) ■ J DY ■ e *Tr(V(z)Y-v(Y))/x {5 21] 



with 

c(Z) = (27rz/A)-* • det V"(Z) ■ A( J^ )} • e ^r(v(z)-v'(z)z)/x (5 22) 

It has been noticed by many authors that the case p = — 1 (i.e. a logarithmic potential 
V(z) = logz) is likely associated with the c = 1 model [ 2q| . We will now proceed to show 
that this is indeed the case. 
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5.3. The Kontsevich-Penner integral 

We have seen that the c = 1 partition function can be succinctly written as a tau- 
function of the Toda Lattice hierarchy 

r(t,t) = (t\S\F) . (5.23) 

For fixed tk we recover a tau-function of the KP hierarchy, which we can study with the 
techniques of the previous subsection. Indeed the operators O n of the minimal models 
should now be compared to the outgoing tachyons of the c = 1 model. 

We want to determine in more detail the element W(t) in the Grassmannian that 
parametrizes this particular orbit of the KP flows. To this end we have to consider the 
state 

oo 

\W(t)) = S- U(t) ■ |0>, U(t) = exp ¥inOi-n ■ (5.24) 

n=l 

We will describe |W(t)) by giving a basis Vk{z]t), k > 0, of one-particle wave-functions. 
First we observe that the operator U(t) acts on the wave-functions z n by simple multipli- 
cation 

U(i) : z n -> exp {^iixi k z~ k ^ ■ z n . (5.25) 
Similarly we have for the action of S a multiplication 

S : z n -> R Pn ■ z n . (5.26) 

We have already seen that the reflection factors R Pn contain all the relevant information 
of the c = 1 matrix model. At radius (3 they can be chosen to be 

n+ k r(i -z^+ 2±i) 

R Pn = {-i»T— \ (1 . £ ■ (5-27) 

1X2-1//) 

(Recall, we are only interested in the perturbative part in (U -1 of this expression.) The 
usual vacuum |0) is spanned by the non-negative powers z k . Therefore the basis elements 
Vk{z;t) of W are simply determined as 



v k (z;t) = c k - S oU(t)z- 
13 



(5.28) 



with a normalization constant Cfc such that Vk(z; 0) = z k . (This corresponds to the normal 
ordering of the ^-matrix in ( |3.9| ).) Since the reflection factor is basically a gamma function, 
the result can be expressed as a Laplace transform 

v k (z;i) = c'{z) -J dyy k - y -W+V-i)/* e H>lv/')> exp QTz^T fc ) (5.29) 

Here the constant c'(z) is given by 

, W = . (5.30, 

These integral representations are of Kontsevich type if and only if (3 = 1, that is, only at 
the self-dual radius. Indeed in that case we have 



(z;t) =c'(z) ■ J dy y k -expifif^y/z -log y + J^tkV k ) (5.31) 



Therefore, following the procedure in ||[|, we can write the following matrix integral 
representation for the generating functional. Define the integral 

a(Z, t) = J dYe ^r { YZ-^V(Y)\ j (5 32) 

where 

V(Y) = -logY + J2hY- k , (5.33) 
and we integrate over positive definite matrices Y . Then we have 

r(t^ = ^, (5-34) 
a{Z, 0) 

with the parameterization 

t n = fi' 1 ■ -TrZ~ n . (5.35) 
n 

Note that with this normalization r(t, 0) = 1, which is appropriate since we consider 

normalized correlation functions. In order to write down the result ( p.34j ) we had to treat 

the incoming and outgoing tachyons very differently, parametrizing the outgoing states 



through (|5.35|) , whereas the coupling coefficients to the incoming states enter the matrix 



integral in a much more straightforward fashion. Equation ( |5.34| ) should be considered 



as an asymptotic expansion in /U -1 , but, for small enough the expansion in these 

variables will be convergent. In some cases, (e.g. the sine-Gordon case considered in |2(J ) 
the expansion has a finite radius of convergence, and as we increase \tk\ beyond the radius 
of convergence we can have phase transitions. 
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5.4- The partition function 

Matrix integrals of the above type have appeared in the work of the mathematicians 
Harer and Zagier [[Z7] and Penner [[20] in their investigations of the Euler characteristic of 
the moduli space -Mg jS of Riemann surfaces with g handles and s punctures. (See [|8| for 



more details on these wonderful calculations.) The double scaling limit of this so-called 
Penner integral was considered by Distler and Vafa [§9| who also speculated on the relation 
with c = 1 string theory. Their work has been followed by a number of papers concerned 
with double scaling limits and multi-critical behaviour of matrix models with logarithmic 
potentials [P3| . All these papers considered essentially the case Z = 1 and i n = 0, in the 
notation of (|5.32| ). 

Distler and Vafa noticed that — after a double scaling limit and an analytic contin- 
uation — the Penner matrix integral could reproduce the c = 1 partition function at the 
self-dual radius (3 = 1. Recall that the free energy at that radius is given by 



^-fief^W^ZV] , (5.36) 



d/j, 2 Jq x \sinha;/2 

and has an expansion 

F = log„ - - log M + g(-l)»^^ 2 - 2s • (5-37) 



(Up to analytic terms in /i.) This makes one wonder whether our result ( p.34j ) can be 
sharpened to give the unnormalized correlation functions. 

To this end let us put the incoming coupling constants ik to zero (and thereby also 
tfc = 0) and take a closer look at the integral 

a(Z) = J dYe ^Tr { YZ^-XozY\ (5 38) 

First of all it has a trivial Z-dependence 

a(Z) = (det Z) N -^ ■ a(l) . (5.39) 

Actually, it is convenient to work with the quantity F defined by 

e f = ( m / / j,)^ 1 e-^ N 2 N / 2 ■ a(l) . (5.40) 
15 



As an asymptotic expansion in \j \x it has the representation 



p J dY ■ e ^ ^k=2 fc 
6 f dY ■ e**™ 2 ' ( ] 



This is known as the Penner integral |2(J and is usually considered in 'Euclidean signature', 
i.e. after analytic continuation \i = iv, v real and positive. 

The quantity F has a beautiful geometrical interpretation, calculating the virtual 
Euler characteristic of moduli space in the open string field theory cell decomposition of 
moduli space. This is essentially the same description of moduli space used by Kontsevich 
0. The expansion of F reads 

oo oo 

^ = EE F Hrf 2 " 2s " S ^- (5-42) 

g=0 8=1 

(Here s > 3 in the case g = 0.) The coefficients are directly related to the Euler numbers 

F SiS = x(M a , s ) . (5.43) 

The 1/fi asymptotics of the integral ( |5.41| ) can be evaluated using the methods de- 
scribed in [|28| to give 

. N-l 

e F = e -W(™)-N/*((- ifl )ii*T(-iiJL)) N {I (l-p/ifx)"-* (5.44) 

^ p=l 

from which one may obtain the formulae: 

'•' 2g(2g - 2 + s) \ s J ' 

It is important that the terms with s = 0, that is, the surfaces without punctures, are 
absent. 

We can explicitly do the summation over s in (|5.42| ) to obtain 

oo 

F = J2 ^ 2 ~ 29 F g (N/ifi) . (5.46) 

9=0 



with 



Fi(i) = -^log(l-i), (5.47) 

~ 3 

F (x) = —5(1 — x) 2 log(l — x) + -x 2 — \x . 
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The double-scaling limit considered by Distler and Vafa in E!J keeps N — i\i fixed, while 



sending N, — > oo (and x — > 1 in (|5.47|) .). This is clearly only possible for imaginary 
H, which is precisely the case they study. However, here we want to consider a simpler 
limit in which \i is kept fixed, but iV tends to infinity. We already mentioned that the 
parameterization ( |5.35| ) only makes sense in this limit. Indeed, the absence of a double 
scaling limit is very much in the spirit of Kontsevich integrals. The contribution for genus 
2 or higher have a smooth limit, as is evident from ( |5.47| ). (Recall, we send x — > oo.) 
However, we have to worry about the genus zero and one pieces, which have to be corrected 
by hand. (This is by the way also true for the double scaling limit.) 

Combining all ingredients we obtain the following final result for the unnormalized 
generating functional for the c = 1 string theory 



r(t,t) =c(Z)- J dYexpi^iTr Y Z~ l - \ogY + J2 tkY ~ 
where the normalization constant is given by 



(5.48) 



c(Z) =e~ lflN (27ti/ij) N / 2 {&etZy»- N 



(5.49) 



The expression ( p.48| ) has a smooth large iV limit. 



6. Other Backgrounds 

The results of the previous sections comprise in principle a calculation of the partition 
function in arbitrary tachyon backgrounds (subject to the equations of motion). The full 
space of classical backgrounds in the theory includes in addition to these excitations of 
the "discrete states" corresponding to global modes like the radius of the ID universe and 
generalizations thereof. Of these, the ones best understood in terms of the matrix model 
are the zero-momentum excitations which are thought to be represented by variations in 
the double-scaled potential. In this section we study the dependence of the amplitudes on 
these extra parameters. We note that in principle the formulation of section three applies 
in arbitrary potentials. What we add here is a study of the variation of the reflection factor 
R q , hence of the partition function, under variations of the potential. 
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6.1. Dependence on (3 

The most obvious parameter is /?, the radius at which we compactify the scalar field 



X. The formulas of section four are valid for arbitrary [3, however as pointed out in [ I7[| , 
correlation functions at different radii are related. The relation is most simply written in 
terms of rescaled couplings t n . Defining 



fi[t n ,t n ;P;tj] = ^F[^-H n ,^-H n -p-^ (6.1) 
so that derivatives of /F yield correlation functions of T q , we have 

J_ _d_ 

F[t n ,i n ;f3;fj] = 2/3 a " N F[^-H{n/P), ^-^(w//?); oo; fj] . (6.2) 

Comparing this with the j3 — > oo limit of the previous calculations is a pretty consis- 
tency check. As an example, set /3 = 1 and consider the two-point function. Computing 
the one-loop graph we find 



dJ-[t n , t n ; 1; fx] -n / • \ fa o\ 
QfQl =V- 2^ R p m R n- Pm =i }^(-iLi-m) n . (6.3) 



m=0 m=0 



Inverting the operator in ( |6.2|) as 



sin( J-) , z-i/2 

— -jp^f= I dsF[n^ fjL + is] (6.4) 



2d{i J -1/2 

we obtain 

(T n T- n ) f3 =oa = i n dx{\-i\i- x) n (6.5) 
./o 

in agreement with the result of [ |10|] . 

6.2. Other zero-momentum modes 

The matrix model naturally suggests candidate representatives of the special states at 
zero X momentum. Operators with the appropriate quantum numbers may be introduced 
as generating variations in the double-scaled potential V{\). Their correlators may thus be 
studied by analysis of the variation of the partition function Z computed above under these 
changes in V. From the definition of I(q, Ai, A2) we can obtain directly constraints on the 
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variation of R q . Essentially these follow upon integration by parts from the linear Gelfand- 
Dikii equation satisfied by a product of Sturm-Liouville eigenfunctions j31|. Explicitly, we 
have 

L qjk R q = k > -1 



= -k(k 2 - 1) -!*— + Aiq{k + 1) A + 2 V Sp (2fc + p + 2) , (6 ' 6) 



where the space of potentials is parametrized by the formal expansion V = J2 n >o s n^ n - 
The operators L q> k for any fixed q are seen to satisfy the commutation relations of (one 
half of) the Virasoro algebra. These were derived from related considerations in [[32]; the 
details of the derivation in the present context as well as the relation to this work appear 
in the appendix. 

Via ( |3.10| ) these imply constraints on the V(A)-dependence of the partition function, 
since this arises only through R q . These however are nontrivial to write down explicitly. 
In particular, we note that they do not seem to fit into the Wi+oo algebra discussed in 



section four. Furthermore, it is easy to see that away from = the identity (^J) ceases 
to hold. Thus perturbing away from the standard background may break the symmetry 
of section four. Further work is required to clarify the relation of the various symmetry 
algebras which appear in this model. 



7. Discussion 

Some remarks on the VF-constraints (|4.10|) are in order. First, it cannot have escaped 
the reader that these constraints are strongly reminiscent of the famous VF-constraints of 
c < 1 models coupled to gravity [§],§[. In these latter models there is only one continuous 
spacetime coordinate and there is only one set of couplings tj rather than the t, i of the 
c = 1 model. Moreover, closer comparison of the identities reveals some important differ- 
ences. For example, at c < 1 the t n couple to an infinite set of gravitational descendents, 
while at c = 1 the t, i couple to gravitational primaries. Nevertheless, a clearer space- 
time interpretation of the c < 1 models will probably emerge from a comparison of these 
identities .! 



5 In [33] proposals for c = 1 flow equations were made by taking the N — > oo limit of the 
Wn constraints of the c < 1 models. It should be noted that, although our equations have some 



similarities to the proposals of [32], they are not equivalent. 
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From the relation of these results to a Kontsevich-type matrix model it appears that 
we have taken a step closer to a unified description of all the c < 1 models along the lines 
proposed by H|§. Moreover, the description ( |5.48| ) of the partition function is a strong 
hint that the c=l correlators have a description in terms of a topological field theory. If 
this is so then the present results provide a direct bridge between a topological field theory 
at the self-dual radius and the local physics of the c = 1 tachyon in the uncompactified 
theory. 

There have been many discussions of Woo symmetry in the c = 1 system. Our con- 
straints are related to the results of []T5| , ^4| -p6| . The other modes of the Woo currents 
appearing in equation ( |4 . 1 0| ) define a set of operators a n {T q ) whose correlation functions 
are determined by the subleading terms proportional to \~\i\~ 2n in the large A asymp- 
totics of the eigenvalue correlators.il These "operators" exist at any radius for X and have 
free fermion representations as fermion bilinears. Their correlators are also given by a 
Toda tau function generalizing that in ( |3.10| ). Note that these operators appear at any 
momentum q and are related to fractional powers of I (or, equivalently, of A). Therefore, 
at generic q they cannot be the special state operators but rather are related to contact 



terms associated to singular geometries created by intersecting macroscopic loops [jloj . At 
integer q the distinction between special states and the a n {T q ) is less clear. We hope to 
return to the subtleties of these contact terms in a future publication. 

The Woo symmetry we have discussed might also be related to the Woo Ward identities 



of [37-43] . In these references the Liouville field is treated as a free field, in other words, 
one works at fi = 0. One should be cautious about identifying these Woo symmetries with 
those of the matrix model. As we have emphasized, the Woo modes of the matrix model 
a n (T q ) are constructed from the tachyon degrees of freedom in distinction to the Woo 



currents of [37-43]. Moreover, our Ward identities are highly nonlinear when expressed 
in terms of the correlation functionsB in contrast to the quadratic identities of [|39|-f£3"|| . 
Finally the ghost sector of the theory is crucial in [|39|-443|], leading to many more "special 



state operators" at given X, cj) momenta than are considered in ]|T5| , |34] -|36|] . Clearly there 
is a certain amount of tension between these two approaches and further work is needed 
to see if these differences are superficial or essential. 



In [ 1 6 1 these operators were denoted cr 2Tl (Og). In [15] it was pointed out that they only have 
contact term interactions. 

7 This is already true at genus zero [18]. 
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We must emphasize that at c = 1 the Woo-constraints are actually somewhat sec- 
ondary, since we have an explicit solution of the appropriate Toda tau function given by 
( p.lO|) . Analogous representations for the c < 1 tau functions (at nontopological points) 
replace the simple operator S by complicated and uncomputable objects like the "star 
operators" of J0| . This is why the Virasoro constraints at c < 1 are essential to the actual 



computation of amplitudes. 

It would be interesting to investigate further the physical properties of these different 
time- dependent backgrounds. In [|18|J26|1 some results along these lines were discussed. Our 



result (|3.10|) should allow a much more complete analysis of the space of time-dependent 
backgrounds in 2D string theory and the various phase transitions occurring as one in- 
creases the coordinates What is needed for further progress is a more effective way 
to compute the tau function (perhaps from the Kontsevich representation) or a deeper 
understanding of the infinite dimensional geometry of the associated Grassmannian. 
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Appendix A. Dependence on the potential 

In this appendix we will derive constraints on the dependence of the free energy upon 
the double-scaled matrix model potential V(X). We restrict attention to variations of the 
potential which preserve the asymptotics V(X) ~ — |A 2 for large A. For convenience, we 
will assume that the asymptotic behavior of V at A — > — oo is such that the eigenfunctions 
decay rapidly enough to render all integrals convergent and all boundary terms negligible. 
To all orders of perturbation theory, of course, these details are irrelevant. 

We wish to compute the change in the correlators ( |3.1| ) under a variation V(A) — ► 
V + 5V(X) preserving the asymptotics. The change in (|3.1| ) is clearly computed from the 
variation of the resolvent I. Thus write 

SI(q, Ai, A 2 ) = (Ai|«5 (— — ] |A 2 ) = - f°° d\5V(X)I(q, X u X)I{q, A, A 2 ) . (A.l) 



H — fx — iq 
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We now recall the calculation of / from []T(| (see appendix A of this work for a detailed 
calculation for particular potentials). For simplicity let q > 0. We will make use of the 
eigenfunctions of H = — V(A) with eigenvalue z = fj, + iq. We have for any V with the 
correct asymptotics two solutions 



a 2 e 



(A.2) 



In terms of these we can write the resolvent quite easily by imposing the boundary condi- 
tions and the defining property (H — z)I(q, X±, X 2 ) = <5(Ai — A2) as in QTOj 



I{q,X 1 ,X 2 ) = -id(X 1 -X 2 ) 



X {z, Ai)x + (z, A 2 ) +R q X (z,Xi)x (z,X 2 ) 



+ (A] 



A 2 ) . 
(A.3) 

The reflection factor R q contains all the effects of the potential, and for the standard V is 
given by (|3.3| ). Inserting this into ( A.l ) and neglecting terms of order SV(Xi j2 ) for large 
Ai, we find that a variation of V yields 



SR, 



dxsvix^iz.xy 



(A.4) 



where %j) = x + + RqX~ is the solution satisfying the boundary conditions at small A. The 
integrand F(z,X) = tfj(z,X) 2 in (|A.4|) satisfies a differential equation |31] following from 
that satisfied by ifi 

F'" — 4(y(A) + z)F' — 2V'F = (A.5) 

where primes denote A differentiation. Let us choose as a convenient set of variations of 
the potential 8V(X) = ee~ iX . Inserting this in ( A.4j) and integrating by parts we findi 



lf-4z£-4£V(-±) + 2V>(-±)}SRq 







(A.6) 



The integration by parts is justified by the limiting conditions we have imposed upon ip 
and 5V. 

Formally expanding V = J2 n>0 s n X n the bounce factor becomes a function of the Sji 
R q = R q [si, s 2 , ■ ■ ■}■ Rewriting 5V as a motion in Sj and inserting the resulting expression 
for SR q in (|A.6| ) we obtain (after shifting sq) 



Lq y kRq 



J q,k 



k > -1 







— — +4i?(Hl)^- + 2V Sp (2Hp + 2)- 

OSk-2 OS k ^ OS p+k 



(A.7) 



8 The similarity of this to the WdW equation of fill is no coincidence; setting z = fi and 
Ai = A2 we find that ( |A.1[ ) is essentially the WdW wavefunction of the cosmological constant. 
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These constraints were obtained in |$2| by different means. We will show that the two 
results are equivalent, but note here that the present derivation has the advantage of 
working with potentials with the correct asymptotics throughout, as well as demonstrating 
explicitly the justification for the various integrations by parts. 



The Virasoro constraints in [52] were obtained as differential equations for the cosmo- 
logical constant one-point function. We have obtained above identical constraints on R q , 
from which one can derive constraints on the partition function. We will now relate the 
two quantities, demonstrating that in fact the two sets of constraints coincide. Begin with 
the formula for the two-point function |l(J[ j: 

(T q T_ q )= f 9 dxR x R* (A.8) 



or its differentiated version: 



§-{T q T. q ) = 2Im[R q R* ] (A.9) 



For q > we may write: 



R q = e *e(/*-Ng;V) (A.10) 

where for E real the pure phase e l ®( E ' ,v ^ is the reflection coefficient for a free fermion of 
energy E in the double-scaled matrix potential V(X). 



We now obtain the specific heat from the limit as q — > 0. As explained in ||4q , [46|j47 
the cosmological constant vertex operator in the c = 1 theory is given by 



(A.ll) 



and is therefore obtained by the limit 



lim -T q = T . (A. 12) 

q^o q 

Taking the limit in ( |A.9| ) we find the leading order begins at q 2 , as expected 
from (|A.12| ). Indeed, quite generally, the low energy theorem of [10 1 shows that for 
n > 2-point functions if k < n momenta qi approach zero the amplitude behaves like 
Ili=i Qi(wp,) k (TYk+i Vq*) m accordance with the expectations of Liouville theory. Thus, 
the apparent - divergence in (|A.12j ) does not appear. From the Liouville point of view this 
may be interpreted as the decoupling of the wrong branch dressing of the vertex operator 
4^,^]. The low energy theorem is more subtle in the case n = 2. In this case the leading 
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order behavior is q + 0(q 2 ). The first term is //-independent and physically sensible, being 
the inverse on-shell propagator at genus zero. The second term of order 0{q 2 ) defines the 
correct zero-momentum two-point function. 

Taking the limit of (|A.8| ) and bearing in mind the above remarks we obtain the 
equation 



(T T ) 



(A.13) 



and hence the nonperturbative one-point function and vacuum energy are 



(To) 



e(n;V)=ilogR(n;V) . 



(A.14) 
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Figure Captions 

Fig. 1. (a) A pictorial version of the integral / for positive momentum, (b) A pictorial 
version of the integral / for negative momentum 

Fig. 2. Incoming and outgoing vertices. The dotted line carrying negative (positive) 
momentum qi should be thought of as an incoming (outgoing) boson with energy 
\qi\. Momentum carried by lines is always conserved as time flows upwards. 

Fig. 3. A real history as a composition of three maps 
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